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Abstract 

The nonlinear equations describing phase ordering dynamics can be closed by as- 
suming the existence of an underlying Gaussian stochastic field which is nonlinearly 
related to the observable order parameter field. We discuss the relation between dif- 
ferent implementations of the Gaussian assumption and consider the limitations of this 
assumption for phase ordering dynamics. The fact that the different approaches gives 
different results is a sign of the breakdown of the Gaussian assumption. We concentrate 
on the non-conserved order parameter case but also touch on the conserved order pa- 
rameter case. We demonstrate that the Gaussian assumption is fundamentally flawed 
in the latter case. 
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I Introduction 



After a system is quenched from the disordered to the ordered phase, domains of the ordered 
phases form and grow. At late stages, it is empirically known that phase ordering process 
obeys dynamical scaling, i.e., the spatial distribution of domains can be described by a single 
time-dependent length, L(t). On this length-scale, the phase ordering process depends only 
on a few general features of the dynamics. Due to the inherently nonlinear nature of the 
dynamics, understanding the phase ordering process remains a challenge Q, g g. Analytic 
progress has been confined to the case of 0(n) component order-parameter in the limit of 
large n [Q, [|, For n < d, where d is the spatial dimension, topological defects become 
important and progress has been limited to dimensional analysis of the defect motion |2|, [7|, [8|] 
and methods by which the nonlinear equations describing the dynamics are "closed" (closure 
approximations) |], [l0|, [ll], 12, 13, 14]. For the scalar order parameter case the simplest 



closure approximation is to assume that the order parameter field <p(r, t) is a Gaussian 
stochastic field ||. This method fails, since, on the scale of the characteristic domain size 
L(t), (j)(r, t) is effectively discontinuous. Instead progress is achieved by assuming that there 



exists an underlying Gaussian field which is nonlinearly related to ip(r,t) fl2| , 13, 14] 



In this paper, we will explore the reliability of closure approximations based on the 
assumption of an underlying Gaussian stochastic field. We first discuss phase ordering 
dynamics without conservation of order parameter. In particular we will discuss the limita- 
tions of this assumption and the relation between different implementations of the Gaussian 
closure. One method of introducing the underlying field is based on the dynamics of ran- 
dom interfaces. We will denote this scheme T. It was first proposed by Ohta, Jasnow 
and Kawasaki (OJK) ]l2j for systems with non-conserved order parameter (NCOP). This 
method was later extended by Ohta and Nozaki [jn| to systems with conserved order pa- 
rameter (COP). 

A second manner of introducing the underlying field is by writing the order parameter 
as ip(r,t) = f(u(r,t)), where / is the interface profile. We will call this method scheme 'B' 
since it relies on the full bulk dynamics. Oono and Puri (OP) |l3[] introduced this scheme 
to remove inconsistencies in the original OJK argument. There are two manners in which 
the bulk closure scheme has been applied. The first (scheme Bii) obtains a closed equation 
for the correlation function of u(r,t). Oono and Puri demonstrated that this provides 
reasonable results for the NCOP case |Tj|. The second (scheme B-0) constructs a closed 
equation for the correlation function of ift(r, t) using the assumed Gaussian nature of the 
underlying field u(r, t). Oono applied the ~Bip closure to the COP case [l(| but found that 
it did not yield the experimentally observed scaling, L(t) ~ i 1//3 . Mazenko has applied the 



Bip closure to both the NCOP [jj and COP cases |17|, |T§. For the NCOP case, Mazenko 
found the B^ closure gives different results from Ohta-Jasnow-Kawasaki and Oono-Puri 
but is also in reasonable agreement with experiment fl4j| . For the COP case, Mazenko found 



that L(t) ~ t 1 / 4 1 17] and included an ad-hoc term to enforce L(t) ~ t 1 / 3 behavior. The 



modified version still contains unphysical features [18] such as the violation of the Tomita 
sum rule [19, pPf] , 
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In this paper we will concentrate on phase ordering dynamics without conservation of 
order parameter but we will also discuss the case in which the order parameter is conserved. 
For the non-conserved order parameter case, the closure schemes reproduce many features 
observed in domain growth, such as dynamical scaling and predictions for the growth of L(t), 
the structure factor IT! 



13, Tj] and the decay of two-time correlations pj], 



. Variations 



of these models have been used to study finite-size effects [24], finite temperature effects 



[13, 25], effects of long-range initial conditions [26], systems with long-range interactions 
|]27|| , and systems with a higher component order parameter |2^, 12^, 3C]. 
predictions of the different closure schemes do not always coincide |31[| . 



However, the 



In Section II, we derive a modified form of the OJK result starting with the interface 
approach (scheme I) correcting for inconsistencies in the original argument [12]. In Section 
III, we discuss the bulk u-closure (scheme Bw) for the NCOP case and demonstrate that 
asymptotically it provides the same result as the interface approach (scheme I). This was 
originally suggested by Oono and Puri using a heuristic argument [O], We next discuss 



the bulk ^-closure (scheme B^). As shown by Liu and Mazenko [31], this closure scheme 
leads to different predictions from the interface approach. We will discuss each approach in 
a parallel manner to emphasize the relation between the approaches. We will show that the 
Bit and Btp closures uses the same approximation so that their different predictions signify 
the breakdown of that approximation. In Section IV we study the approximate Gaussian 
nature of u(r, t) through numerical simulations. We find that the single point probability 
distribution function P(u) decays as a Gaussian at the tails but decays slower than Gaussian 
near u = 0. In Section V we present our findings for the conserved order parameter case. 
We find that, in this case, the Gaussian closure is fundamentally flawed. In Section VI we 
summarize our findings. 



II Interface Approach 

In this section, we derive a modified form of the Ohta-Jasnow-Kawasaki and Oono-Puri 
results demonstrating the features of the interface scheme. We assume the dynamics are 
described by the time-dependent Ginzburg-Landau (TDGL) equation 

^ = -WW) + |vV, (2-1) 

where £ is the interfacial width, ip(r,t) is the scalar order parameter and iib{'4 ) ) is the 
portion of the local chemical potential which contains no gradient terms. We assume [xb 
is an odd function of ift, the equilibrium values of ij; are ±1; and there exists a solution, 
tp = f(z), corresponding to a stationary planar interface at z = 0, i.e., f(z) obeys 

d 2 f 

o = -M/(*)) + V^- 

For example, if lib = —ip + tp 3 then f(z) = tanh(z/£). However, other than the requirement 
that f(z) increases monotonically from —1 to 1 over a length-scale ^, the exact form f(z) is 
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unimportant [jl3| , p^] . This universality is closely related to the universality in the dispersion 
relation around the interface. 



At late times, the width £ is small relative to the characteristic domain size L(t), and 
the domain growth is determined by the motion of the sharp interfaces. The interface 
dynamics can be derived from the TDGL @, |, @, 

e 2 

u n (r,i) = -^-K(r,t). (2.2) 

Here the normal velocity of the interface, u n , is defined as positive when the 'minus' phase 
moves into the 'plus' phase, the normal n points into the 'plus' phase; and the local curva- 
ture, k = V • n, is positive for a bump of the 'minus' phase into the 'plus' phase. 

The essence of the interface method is to rewrite the interfacial equation in terms of an 
indicator field u(r, t). The indicator field is defined so that u > (u < 0) in the plus (minus) 
phase, u = at the interface, and, near the interface, \u\ is the distance to the interface, so 



that, Vu(r, t) = n(r, t) near u = 12, 13]. The main motivation for introducing u(r,t) is 
that it remain continuous on all length-scales, while ip(r, t) is effectively discontinuous upon 
rescaling distances by L(t). Therefore simple decoupling approximations should be more 
trustworthy for the indicator field u than for order parameter field tp. 

In terms of u the interfacial dynamics (Eq. (|2.2| )) is 

« = y v v (2 - 3 > 

However, Eq. ( |2.3| ) can only hold at u = 0, since, if it were true in the bulk, the condition 
that |Vit| = 1 near the interface would be violated ]l~3[ . We assume that the |Vu| condition 
can be met by extending Eq. (|2.3| ) into the bulk by adding a "Lagrange multiplier function" 
P(u,Vu) = P(u,Vu)u to the RHS of Eq. (O), 



^ = f [v 2 u + P(u,Vu)u . (2.1) 



The function P has the following properties: P(0, V«) = 0, since Eq. (2.3) must be recovered 
at u = 0, Due to the symmetry of the TDGL, P(u, Vit) is an odd function of u and an 
isotropic function of derivatives of u. Since the interface dynamics depends only on local 
properties of the interface, P(m, Vm) is assumed to be local, i.e., P(u,Vu) depends only 
on a finite number of derivatives of u. Using heuristic arguments, Oono and Puri replaced 
P(u, Vu) by a function of time chosen to maintain the equilibrium interface thickness [|1J. In 
the following section, we demonstrate that P(u, Vit) can be obtained exactly from the TDGL 
using the bulk u-closure (Bn approach) and therefore schemes I and Bu are asymptotically 
equivalent for the non-conserved order parameter case. 



d{mu 2 ) 



From Eq. ([2.4]), the two-point correlation function, (uiu^), obeys 

(V?(uiu2> + {P{u h Vi«i)n 2 )) , (2.5) 



dh 2 
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where Ui = u(r.j,tj) and t\ ^ t 2 . Note that the local constraint |Vit| = 1 at the interface 
forces ( u 2 ) to grow as L 2 , where Li = L{tj). 

We now need to make an assumption concerning the statistics of {u(r, t)}. The simplest 
assumption is that u is a Gaussian stochastic field. As noted previously, the rationale is 
that u is 'equicontinuous' for all length scales, so simple decoupling scheme may be more 
trustworthy for it than for the original order parameter field. 

For any Euclidean symmetric Gaussian stochastic field, {P{u\, Viti) u\u 2 ) has a simple 

form 

(P(m, Vm) uiu 2 ) = p({ u\ ), V 2 ( u\U\ )) { u%u 2 ), (2.6) 

where V 2 {u\Ui ) = rim r /_ >ri V 2 (it(ri, ti)u(r' , ti)). Note that the function p does not depend 
on t 2 or ri — v 2 so we can replace p as a function of t\: 

(P(ui, Viti) uiu 2 ) = p(h) ( u\u 2 ). (2.7) 

For t\ / t 2 Eq. (|2.5|) becomes 



d(uiu 2 ) _ C 
dti ~ 2 

For t\ = t 2 = t, we have 

dg(r, t) 



2 . . 

r (Vf(«i«2)+P(tl)(«l«2>j • (2i 
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^(V 2 9 (^i)+P(i)5(r,t)) , (2.£ 



where g(r,t) = (u(r,t)u(0,t)). The Gaussian approximation means that the detailed local 
constraint \Vu\ = 1 is no longer met but is replaced by a global constraint (u 2 ) ~ I? . This 
requires that i 2 p{t) = (d + 2)/(2t), where d is the spatial dimension. This is exactly the 
form of p{t) chosen in so that the interface thickness is time independent. 



The final step is to obtain the correlation function of ip{r,t) from u(r,t). Near the 
interface ifi will be close to the planar interfacial profile so that ip(r,t) = f(u(r,t)). At long 
times, the exact form of / should be irrelevant, (as long as / approaches sgn(u) as £ — > 0) 
]. Oono and Puri chose a form that simplifies the Gaussian integrals, 



1 ( u 2 



f(u) = sgnH * {2<2)1/2 exp y-^j ( 2 - 10 ) 

where * is the convolution. They find 

2 / ( u\u 2 ) \ 

f 7 12W = -arc 8 m^ ((t(;)+{2)1/2({t( , )+al/2 j, (2.11) 

where C\ 2 {r) = {ipiip 2 ) and r = \r% — r 2 \. For any reasonable form of f(u), a WKB-type 
argument can be used to show that, C\ 2 {r) approaches Eq. (|2.11|) when (u 2 ) 3> £ 2 . 

In the original OJK analysis (as well as OP) there is a factor (d — l)/d in front of 
the diffusive term in the evolution equation for u. This difference is due to the lack of 
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local constraint, |Vu| = 1, in OJK. (For OP, this is due to both its eclectic and logically 
opaque nature and the lack of this exact constraint.) With this exception, Eq. (|2.9| ) with 
p(t) = (d + 2)/(2£ 2 i) together with Eq. fl2Tl| ) is the main result in Ref. @. The result 



of OJK can be obtained from this approach asymptotically without the difficulties in the 
original argument. In the original OJK derivation u is not self-consistently defined, since / 
was chosen to be sgn(ti) and |Vu| = 1 was not even approximately enforced [12, 



The interface approach predicts the following |l2|, (1) dynamical scaling with 

L ~ t 1 / 2 , (2) the equal-time correlation function, C(r,t) decays as a Gaussian at large r, (3) 
the scattering function, Sk, displays Porod's law, Sk ~ p5|, 19], for the wavevectors 



k in the range L 1 <C k <C ^ 1 and (4) Tomita's sum rule is obeyed [19|, i.e., in the scaling 
limit, d m C/dx rn \ x=r /£ j ^Q+ = for any even m. For the two-time behavior, the interface 
approach predicts C\2(0) ~ (t\jt<2)~ d l 2 for £i / £2 <C 1 and the autocorrelation function of the 
Fourier transforms of ip decays as a stretched exponential, ln( V'k^iJV'-k^) ) ~ — (^/ti) 1 ^ 2 , 
for L 2 » Li [pll. 



An alternative closure method is the one due to Kawasaki, Yalabik and Gunton fll lfl . 
This method tries to sum a diverging series whose n-th term is of order exp(n7t) with 
jt 3> 1. Although the result is finite and the scaling results are the same as that of 
Eq. (|2.S|), this diverging behavior causes the ratio ^/L{t) to vanish exponentially fast, an 
absurdity caused perhaps by the strongly diverging nature of the series. 



Ill Bulk Approach 

III. A Bulk w-closure (Scheme Bm) 

In this subsection, we discuss the bulk u-closure (Bn). We demonstrate that the results 
of the interface approach can be obtained without the intermediate step of the interface 
description. 

We introduce an auxiliary field u(r, t) by 

VM) =/(u(r,t)). (3.1) 
The planar interface solution, /, obeys 

^ (/) = -2d^' 

At late stages, the interfacial profiles will be very close to the planar profile. Therefore this 
nonlinear mapping is exactly of the form required to enforce u = at the interface and 
V« = n near the interface. 



In terms of f(u), the TDGL is 

(3.2) 



at 2 



^4-v 2 / 

du 2 
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Applying the chain rule for df /dt and V 2 /, we obtain an expression for du/dt, 



du 



-.2 

'-[V 2 u+(l-\Vu\ 2 )Q(u) 



(3.3) 



where Q(u) = — (df/du)~ 1 d 2 f/du 2 . For \xb = — V' + V' 3 we nn d f{u) = tanh(ii/£) and 
Q(u) = 2f(u). In general, Q(0) = 0, Q(u) is an odd function of u and Q{u) approaches 
d/j,B/dip\^ eq (a finite constant) as u — > oo. The result is that Q(u) must be proportional to 
sgn(u) in the limit £ — > 0. 



Equation (^j) is exactly the form required to extend the interface equation into the 
bulk using the interface approach (Eq. (2.4)) of the previous section. We can now identify 
the previously unknown 'Lagrange multiplier' function, P(n, Vn) in Eq. (2^), with (1 — 
\X7u\ 2 )Q(u). We repeat the steps of the interface approach with the difference that we now 
have an explicit expression for P(u, Vit) 



(3.4) 



From Eq. ( p.4j ) we obtain the expression for the correlation function {u\U2j 
d(uiu 2 ) 



Ot 



r v>iU2> + <(i-ivi«ir)Q(ui)u2) 



where we restrict the discussion to t\ = t% = t. The more general case is easily obtained. 

To this point there are no approximations. To proceed further we again assume that 
{u(r,t)} is a Gaussian stochastic field. Just as in scheme I, this assumption is totally 
uncontrolled. In the same spirit as before, we assume that explicit forms of / and Q are 
unimportant as long as they approach sgn(u) and 2sgn(n) in the limit of £ — > 0. Therefore 
we choose a form of Q which simplifies the Gaussian integrals, 



i r+oo i 
Q{u) = — I dw— exp 

7T J-oo W 



2, ,2 



UjJU 



That is, Q(u) is the sign function mollified by a Gaussian function. Equation (| 
be computed explicitly as 



dg(r,t) 
dt 



V 2 g(r,t) + 



1 + V 2 g(0,t) 
Vvr(e+5(0,t))/2 



(3.5) 



can now 



(3.6) 



where, as in the previous section, g(r, t) = ( uiu 2 ) for t\ = t 2 
equation exactly the same form as Eq. (^9|) 



dg(r, t) 



dt 



where 



Pit) 



V 2 g{r,t)+p(t)g(r,t) 

1 + V 2 g(0,t) 
V^(e 2 + 5(0,t))/2 



t. That is, we have an 
(3.7) 
(3.8) 



is now an explicitly known function of t. 
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We are left with a closed equation for g(r, t) which can be solved as an initial value 
problem. However, it is more instructive to consider the scaling limit. We define the charac- 
teristic length-scale by L 2 = g(0, t)/(-d 2 g(r, t)/dr 2 \ r ^ ) so that V 2 g(0,t) = -dg(0, t)/L 2 . 
With this definition we can show that L 2 2^H. Substituting this in Eq. (U) we find 



dg(0,t) _ d _ fn ^ t 2 



dt 2t 



g(0,t) + ?p(t)g(0,t), (3.9) 



with p now given by 



e-dg(0,t)/(2t) 
VT(e + 5(0,t))/2- 

The only meaningful asymptotic solution of this equation is g(0,t) ~ 2£ t 2 t/d + 0(1). This 
implies that asymptotically £ 2 p{t) = (d + 2)/{2t), i.e., the same p(t) as that given by the 
interface approach. We also find that this is the asymptotic form if we solve the initial value 
problem using the explicit expression (Eq. (|3.8| )) for p(t). 



Therefore we have demonstrated that the results of the interface approach can be 
obtained directly from the time-dependent Ginzburg-Landau equation without explicitly 
using the interface description. In particular we have shown that the previously unknown 
functions, whose properties were obtained by scaling arguments and intuitive physical re- 
quirements, can be obtained more directly and, furthermore, is exactly the form required 
by the scaling arguments. 



III.B Bulk ^-closure (Scheme B-0) 

In the bulk ^-closure the auxiliary field u(r, t) is introduced in the same way as in the bulk 
it-closure, i.e., ip(r,t) = f(u(r,t)), where / is the planar interfacial profile. A dynamical 
equation is then obtained directly for ( il>(ri,ti) VK 1 ^?^) ) under the assumption that 
{u(r,t)} is a Gaussian random field. Mazenko applied this scheme to the non-conserved 
order parameter case [14|. We summarize the argument below. 



The TDGL equation in terms of ip = f(u) is given by Eq. ( |3.2| ), 

V 2 / ) . (3.11) 



df efd 2 f _ 2 . 



dt 2 \du 2 

The expression for the correlation function ( ipi ip2 ) = ( fl /2 ) becomes 

9(/iA> _*(i<Pf 



m f 2 ( ( ^r /2> - v2</l/2) ]' (3 ' 12> 

where fi = /(n(rj,tj)). Here, we again let t\ = ti = t for simplicity. An analogous equation 
can be written for t\ ^ £2- 

To this point there has been no approximations. To proceed further, u is again assume 
to be a Gaussian stochastic field. Using the properties of Gaussian variables, one finds 
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( (d 2 fi/du\)u2 ) = d{ /1/2 )/d(u 2 ). This gives a closed equation for C(r,t) = (i/j(r,t)i/j(0,t)). 
The equation can then be solved numerically as an initial value problem. For later times, 
the exact form of / is irrelevant and we can use the / given in Eq. ( 2.10| ). Performing the 



Gaussian integrals, we find, dC/d{u 2 ) = — 2tan(7rC/2)/2/(7r((n 2 ) +£ )), and 

i*7(r,*)_ 2 +f!))t j^ +v!(;(rjf) . 



£ 2 ^ tt((u 2 

For (u 2 ) S> ^ 2 a scaling form can be assumed. A nonlinear eigenvalue problem is then found 
which can be solved numerically to obtain the correlation function C(r,t) Wl 



This closure yields |L4|]: (1) dynamical scaling with L ~ i 1 / 2 , (2) that the quasi-elastic 
scattering intensity obeys Porod's law at large k, (3) that Tomita's sum rule is obeyed fl9f , 
and (4) that the real space correlation function decays slightly slower than a Gaussian with 
C(x) ~ x a exp(— bx 2 ) where a > depends on A and, hence, d. For the two-time behavior, 
Liu and Mazenko predicts that 6*12(0) ~ {L\/L2) h for L\ <C L2 where b ~ 1.28 for d = 2 



and b = 1.63 for d = 3 [23]. This closure also predicts that the auto-correlation function for 



the Fourier transforms will decay exponentially in £2/^1 for ^2 3> t\ [23]. 



III.C The relation between the different approaches 

From the above discussion the bulk u-closure and interface approaches are equivalent. We 
also see that the Bu-closure and the B-0-closure uses exactly the same assumption: tp = f(u) 
with u assumed to be Gaussian. However, we also see that there are discrepancies in the 
predictions of the two approaches. These discrepancies can be regarded as an indicator of 



the limit of the reliability of the Gaussian assumption. Experimentally [36] there is general 
consensus with regards to the features in which both approaches agree, i.e., L ~ i 1 / 2 , 
Porod's law [35], and Tomita's sum rule [19], although these are kinematic consequences. 

To investigate the discrepancies between the different models, one must rely on numer- 
ical studies, either kinetic Ising model simulations with Glauber dynamics or cell dynamical 



schemes (CDS) corresponding to the time-dependent Ginzburg-Landau equation |52|, 37fl . 
Humayun and Bray pGfl , and Shinozaki [|38| find that the equal time correlation function 
C(x) obtained by a CDS is better fitted by the OJK result, although, C(x) decays faster 
than both predictions. Liu and Mazenko studied the autocorrelation function 6*12(0) and 
found that the decay exponent b is approximately 1.25 for d = 2 and for d = 3, preliminary 
results indicate b ~ 1.8, i,e., the decay is faster than the OJK prediction of d/2 and more 
in agreement with the prediction of the B^ closure. On the other hand for smaller ratios of 
times £2/^1 < 10, the OJK result seems to give a very good zero parameter fit to 6*12(0) pl| . 
The numerical data for the autocorrelation of the Fourier transform ipk 3£] has not yet been 



sufficient to compare the results of the two closures. In this case the bulk n-closure pre- 
dicts that {ipk(ti)ip-k(t2)} decays as a stretched exponential for £2 ^ t\ [^] while the bulk 
^-closure predicts an exponential decay p3[| . Nevertheless for many features, the Gaussian 
closures give reasonable agreement with experiment. A broad statement would be that, 
while B^ gives better fits to the two-time behavior, schemes Bn and, hence, the scheme I, 
gives better fits to the equal time behavior. 
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The relation between the closures is further clarified in the d — ► oo limit. It was 



shown by Liu and Mazenko [31] that the scaling functions obtained from OJK and the 
Bip closure coincide in the limit of d — > oo. Here we show more explicitly that the two 
closures are equivalent in this limit. In the limit of t —* oo {ipifo) depends only on the ratio 
a(r,t) = g(r,t)/g(0,t) for t = t\ = t%. The evolution equation for a(r,t) in the Bu closure 
can be obtained from Eq. Q3.6| ). In the limit of large t this becomes 

For the ~Bip closure we can obtain the dynamical expression for a(r,t) from Eq. ( |3.13| ), 

da(r,t) _ e (^L +v 2 a(r,t)+«(r,t) |Va(r,t))| 2 V (3.15) 



dt s V {u 2 ) 

We can examine the d — ► oo limit by rescaling lengths by x = r/yfd. Eq. ( [3.14[ ) becomes 

d<x{x,t) = 2 ( a{x,t) | l d 2 a(x,t) | d-lda(x,t) \ 
dt \ (u 2 ) d dx 2 d x dx J 

While Eq. ( |3.15| ) becomes 



da(x,t) 2 ( cn(x,i) ld 2 a(x,t) d—lda(x,t) a(x,t)fd x 



a 



2^ 



dt * \ (u 2 ) + d dx 2 + dx dx + d \ dt ) J ' ^ 3 ' 17 ^ 

Taking the limit of d — ► oo, we find that the equations for a(r, t) become identical. Although 
this is not a proof, it does indicate that the Gaussian approximation may be valid in the 
limit of d — > oo. 



IV Numerical Results 



In the previous discussion we have shown that the {u(r, t)} is not a Gaussian random field. 
However, we have also argued that this assumption may be a reasonable first approximation. 
In this section we study the statistics of {u(r,t)} directly through a numerical updating of 
Eq. (|3.6| ), As noted previously, no approximations are needed to proceed from the TDGL 
(Eq. (|2,1|)) to Eq. (|3.3|) so this is equivalent to a simulation of the TDGL equation. 



For numerical efficiency we choose Q(u) = 2 if u > 1, Q(u) = 2uifl>u> —1 and 
Q(u) = —2 if u < — 1. This effectively approaches 2 sgn(-u) in the limit of L(t) 3> 1. We 
discretize the system with mesh size 8x = 1.0 and time steps 5t = .05. To reduce lattice 



effects we used a sphericalized Laplacian as described in Ref. [32]. For these large time 
and space steps the update corresponds to a CDS [p^] . The simulation was performed on 
800 x 800 lattices with periodic boundary conditions, and repeated on 400 x 400 lattices to 
check for finite size effects. The results for n = 400 and n = 800 begins to deviate at about 
t = 400 indicating that the data for n = 800 is not affected by size problems. An average 
was taken over 18 independent initial conditions for the larger lattice. 
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Figure 1 shows a plot of L(t) 2 vs. (u(r,t) 2 ) for t = 25,50,100,200 and 400. For 
L(t) we use the inverse interfacial density. The line has a slope of unity. We see that 
Lit) 2 ~ (u(r,t) 2 } in agreement with the arguments above. We next calculate the single 
point probability function P(u,t) for each time. Figure 2 shows In P(u 2 / (u 2 )) vs. u 2 /(u 2 } 
for t = 50, 100, 200 and 400. It is clear that the probability distribution scales during this 
time range. We also observe that the tail of the probability distribution decays in a Gaussian 
manner but the region near u = flatter than that for the Gaussian distribution. Figure 
3 shows the same data plotted with — In (— In P(u 2 /(it 2 ))) vs. \nu 2 / (u 2 ) . The straight line 
has a slope of 2. Figure 4 shows the flatness {u }/({u 2 } 2 ). The times are the same as that 
of Figure 1 plus the point at t = 0. The flatness is 3 at t = since the initial distribution 
is Gaussian. For larger times we find that flatness is somewhat less than that expected for 
a Gaussian distribution. 

From our numerical result we find that the single point probability distribution of u 
is approximately Gaussian at the tails. However, the deviation from Gaussian behavior 
near it = is very important since the location of the interface is at u = which we have 
assumed controls the dynamics. We also note that we have only looked at the single point 
distribution function. To test the full Gaussianness of {u(r, t)} we must look also at the 
two-point and two-time distributions. 



V Conserved order parameter case 

The bulk ip closure (B^-scheme) approach was originally introduced to study systems with 
conserved order parameter (spinodal decomposition). Here we summarize our study of the 
conserved order parameter case emphasizing the difficulties of this approach. 

The first important difficulty is that one cannot simply use the nonlinear mapping 
ip(r, t) = f(u(r, t)) to define the indicator field. This is because the condition that \ip(r, t)\ < 
ipeq is not met so that the mapping is not invertible. This is easily seen since, for the 
conserved case, local equilibrium near the interface requires that the deviation of ip from 
its planar interface value is proportional to the local curvature. However an indicator field 
can still be introduced using 

V>(r,t) = /(u(r,t)) + 0(r,t) 

where the <j) field accounts for the deviation from the planar interfacial profile and the u 
field has the same properties as in the nonconserved case. One can then make the Gaussian 
assumption for the u field with the <f> field coupled in such a way as to enforce the L(t) ~ i 1 / 3 



growth [ 18 1 . 



However, even with this extension we find that the Gaussian assumption is incompatible 
with conserved order parameter dynamics. Since the extra field (ft is of order 1/L we can 
neglect its direct effects in the correlation function in the scaling limit. We note that Eq. 
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2.11), the relation between (iptp) and (uu) 



ni +\ 2 ■ ( ( u i n 2 
C (r, t) = — arcsm 



7T \ (u 2 ) 

is independent of the dynamics and is true as long as u(r, t) is a Gaussian stochastic field. 
Therefore, assuming u is a Gaussian field, we can invert this relation to obtain ( u\U2 )/{ u 2 ) 
from an empirically obtained C(r,t). Figure 5 shows the spectral density (ufc(i)u_fc(t)) 
obtained in this manner using C(r, t) from a very accurate three dimensional CDS simulation 
of spinodal decomposition |40| ]. We find that the spectral density {uk(t)u-k(i)} becomes 
significantly negative at small wavenumbers q = kL(t) < 0.5. (The peak of {ipk(t)i>-k(t)) 
occurs at approximately q = 1.) Since the spectral density must be positive definite, we can 
say that the Gaussian assumption is not a reasonable description of the large length-scale 
behavior. In addition we find that (inset of Figure 5) there is also a violation at positivity 
in the very important range of q from approximately 1.5 to 4. This corresponds to the 
structure at wavenumbers just above that of the peak of the scattering intensity. We can 
conclude that the bulk closure approaches are inherently flawed for the conserved order 
parameter case. 

Since interface approaches also use the Gaussian assumption and the same relation 



Eq.( 2.11 ), they are also flawed. However, Ohta and Nozaki's modest success justifies a 
more systematic study of the interface approach. This is especially the case since the Ohta- 
Nozaki approach fixes the growth exponent to be 1/3 by an intuitive but rather ad hoc 



manner |15|| . Here we briefly summarize our interface attempt and its limitation. 

Our starting point is the interface dynamic equation in terms of the u field 

G * 5(u)d t u = V 2 u, (5.1) 

where G is the Laplacian Green's function, i.e., G = — V~ 2 E|. As in the non-conserved 
case, this equation is correct only at the interface. We assume the following bulk extension 
of the interface equation: 

G * 6(u)d t u = JV 2 u + Q, (5.2) 

where J is a function of u such that J = 1 at the interface, and Q is a functional of u and its 
derivatives up to the second order with Q = at the interface. The extension is motivated 



by an analogous idea behind Eq. ( |2.4j ). Applying the Laplacian to Eq. (5.2), and assuming 
that u is a Gaussian stochastic field, we find (after some algebra) 

\ To \n + A ^r,t) = -R(t)V*g(r,t)-V 2 Q(t)g(r,t)+P(t), (5.3) 
2 y/2TTg(0,t) 

where g(r,t) = (u(r,i)u(0,i)} and P, Q and R are yet unspecified functions of time only. 
We now assume that there is a scaling regime and determine the forms of R(t), Q(t) and 
P(t) necessary for a scaling solution to exist. From the definition of u we require (u 2 ) ~ L 2 
for large times, so that, in the scaling limit, g(r, t) must be of the form 

g(r,t) =L(t) 2 g(r/L(t)), (5.4) 
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where g(0) = 1. Rescaling x = r/L(t), we get Eq. (5.3) with Eq. ( p.4| ) as 



X ^2g(x)-x^j = -R(t)V 4 g(x)-L(t) 2 V 2 Q(t)g(x)+L(t) 4 P(t)g(x), (5.5) 

where A = L'(t)L(t) 2 /2\/2~k which must be time-independent asymptotically. The coeffi- 
cients R,L 2 Q and L 4 P must converge to a non-zero constant in the t —* oo limit; if they 
diverge, we get physically absurd results, while the same is true if these coefficients vanish. 
Hence, asymptotically Eq. (|5.5[) becomes 



X(2g(x) - x^l) + V 4 g(x) + AV 2 g(x) - Bg{x) = 0, (5.6) 

where we have rescaled A to get rid of the numerical coefficient in front of the double 
Laplacian, and A and B are constant. For g(0) to be finite, B = (2 + d)X is required. 
Hence, in the /c-space, we arrive at 

A^| = (-fc 3 + %, (5.7) 

where A is scaled out, which must be positive and A remains an unspecified constant. The 
resultant equation is very similar to the Ohta-Nozaki equation, and has the same defect, 
although in our case L ~ t 1 / 3 follows from our starting point of the interface equation. 

Given these caveats we fit the empirically obtained C(r,t), using the g(r,t) obeying 
Eq. (|5.7|). We find that A = 0.013 gives the best fit to C(r). Figure 6 shows that the fit is 
very good up to approximately the second zero of C(r). This is further shown in the inset 
in which (r / L) 2 C(r) is plotted to show that goodness of the fit is not simply because C(r) 
is small. As noted previously, the Gaussian assumption is invalid for long length scales and 
the fit becomes increasingly worse for larger x. The same holds for the scattering intensity. 
We find a very good fit for q > 0.5 but the conservation law is violated due to the invalidity 
of the Gaussian field assumption at small k. 

To conclude, we find that the Gaussian assumption can give a nontrivial fit to a limited 
range of length-scales for the correlation function and wavenumbers for the scattering inten- 
sity. However, given the fundamental flaw in the Gaussian assumption at larger distances, 
we feel that there is no point in making a more concerted effort f42j|. 



VI Discussion and Summary 

We have discussed the relation between different closure approximations for phase ordering 
without conservation of order parameter. These closure approximations are based on the 
assumption of an underlying Gaussian stochastic field u(r, t). We discuss two general meth- 
ods, the interface approach in which a dynamical equation for u(r, t) is obtained from the 
interfacial dynamics and the bulk approaches in which the u(r, t) is defined by the nonlinear 
relation ip(r,t) = f(u(r,t)), where / is the planar interfacial profile. The bulk approach is 
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further subdivided into the bulk w-closure and bulk ^-closure in which dynamical equations 
are obtained for the correlation function of u and tp, respectively. E| 



We have derived a modified form of the original interface approach [12] and shown 
that it is completely equivalent to the bulk u-closure. We then have demonstrated that the 
only assumptions of the bulk -^-closure and bulk -u-closure is that {u(r, t)} is a Gaussian 
stochastic field. Our conclusion is that the discrepancies in the predictions of the bulk 
■u-closure and bulk -^-closure is due to a breakdown of that approximation. 

We have also discussed the Gaussian closure for the conserved case. We have shown 
that the Gaussian approximation is more fundamentally flawed in this case. However, the 
interface closure approach still leads to a nontrivial fit of the correlation function C(r,t) 
and scattering intensities Sk(t) except at large scaled distances r/L(t) or at small scaled 
wavenumbers kL(t). 
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Figure 1: A log-log plot of (u 2 ) vs. L(t) 2 from the 800 x 800 simulations. The solid line is a 
slope of 1. The statistical uncertainties are smaller than the symbol sizes. L(t) is the inverse 
interfacial density. We find that (u 2 ) grows as L 2 as predicted by the interface approach. 
A fit to the form (u 2 ) = aL 2b + c gives b = 1.03. 



Figure 2: A plot of the single point probability distribution In P(u 2 /(u 2 ) ,t) versus u 2 /(u 2 } 
for t = 50, 100, 200 and 400. Representative error bars are shown. This plot indicates 
that the tail of the distribution function decays as a Gaussian but there is a regime for 
u 2 /(u 2 } < .5 which decays slower than predicted by the tails. 



Figure 3: The same data as in Fig. 2 plotted in the form — ln(— lnP(u 2 / (u 2 ) , t)) vs. 
ln(u 2 /(u 2 )). The line has a slope of —1 indicating the Gaussian nature of the tail. 



Figure 4: The flatness (u 4 ) / (u 2 ) 2 versus time. The distribution for the initial condition is 
Gaussian so that at t = the flatness is 3. For larger times the flatness is somewhat smaller 
than that of a Gaussian distribution. 



Figure 5: Plot of the spectral density (uk(i)u-k(t)) from the (ipk(t)' l P-k('t)} obtained under 
the Gaussian assumption from the 3-d spinodal decomposition simulation (dashed line). 
The inset is a blow up of the spectral density for values of q just above the peak. We find 
that there is a violation of positivity both q < .5 and 1.5 < q < 4, indicating that the 
Gaussian assumption is clearly invalid for the conserved order parameter case. 



Figure 6: Plot of the real space correlation function C(r, t) from the 3-d CDS simulations 
(dash line) versus the result using the Gaussian assumption with the Bip scheme (solid line) . 
The inset shows (r /L) 2 C(r,t) for the same range of r/L. We find very good fits up to the 
second zero in the correlation function. 
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